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Abstract 

The H dihyperon (DH) is studied in the framework of the SU(3) chiral quark model. It 
is shown that except the a chiral field, the overall effect of the other SU(3) chiral fields is 
destructive in forming a stable DH. The resultant mass of DH in a three coupled channel 
calculation is ranged from 2225 MeV to 2234 MeV. 
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In 1977, Jaffe predicted DH fiL a six-quark state with strangeness (s) being -2 and 
J p = + (S = 0, T = 0), by using a simple color magnetic interaction in the MIT bag model. 
Since then, many theoretical [Q, [|, [3], ^, ||, ffl] and experimental [|| (and references therein) 
efforts have been devoted to the DH study. There were so many theoretical predictions of 
DH which are quite different in different models. For instance, by using the MIT bag model, 
Jaffe gave a binding energy of about 80MeV below the two A threshold (E\\) |IJ; in terms of 
the Skyrme model, Balachandran et al. even showed a larger binding energy of about several 
hundred MeV ||]; by employing the cluster model, Yazaki et al. predicted the energy of DH 



from about WMeV above E\\ to about lOMeV below £"aa |§ |14|| ; also, by using the cluster 
model but different interaction with Yazaki's, Straub et al. announced a binding energy 
around 20MeV J|]; in terms of Quantum Chromodynamics((5C-D) sum rules, Kodama et al. 
gave a binding energy around AQMeV although the error bar was quite large; considering 
mutually the two-cluster and six-quark cluster configurations, Wolfe et al. obtained a deeply 
bound DH with a binding energy about several hundred MeV Q; and by employing a quark 
model without the one-gluon-exchange interaction (OGE), Glozman et al. announced the 
non-existence of a bound DH On the other hand, there is no experimental evidence 
showing the existence of DH up to now. The provided lower limit of the DH mass is about 
2200MeU g [ilj. 

The reason for carrying out such researches is straightforward. According to the feature of 
the color magnetic interaction (CMI) in the one-gluon-exchange potential (OGE), when the 
strangeness of the system concerned is equal to -2, S=0 and (\^)f = (0 0), the expectation 
value of CMI presents more attractive feature than those contributed by two A baryons. 
As a consequence, six quarks could be squeezed in a small region, the typical short-range 
QCD behavior would be demonstrated and some new physics might be revealed. Therefore, 
studying DH will be rather significant in understanding the quark characteristics of the wave 
function of the multiquark system and the short-range behavior of the QCD theory. 

On the other hand, till now, there still exist some uncertainties in the nucleon-hyperon 
(NY) interaction on the baryon-meson degrees of freedom, especially in the short-range part, 
so as the hyperon-hyperon (YY) interaction and the prediction of DH. Therefore, the next 
generation, the quark-gluon degrees of freedom, may be a more effective base to establish the 
relations among the nucleon-nucleon (NN), NY and YY interactions, and consequently to 
give a more reliable prediction of DH. As well known, most nuclear phenomena are just of the 
low energy approximation of QCD. There exist lots of nonperturbative QCD (NPQCD) 
effects. Unfortunately, nowadays one still cannot solve NPQCD directly, and has to employs 
certain QCD inspired models. The SU(3) chiral quark model is just one of the most successful 
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ones. In that model, the couplings between chiral fields and quark fields were introduced to 
describe the short- and medium-range NPQCD effects, and a more reasonable quark-quark 
interaction V q - q (q = u,d, and was obtained. With that V q - q , one could mutually 

describe the experimental NN scattering phase shifts, available YN scattering cross sections 



1 11] and some properties of single baryons including the empirical masses of single baryon 



ground states [12, 13]. Extrapolating that model to the s = —2 system, one can study the 
YY interaction, double strangeness hypernuclei, DH and etc.. In this letter, we would choose 
DH as a target, because it is a simplest case with two strangeness and its structure is relative 
simple so that the interaction can be preliminarily examined by the present experimental 
finding although it is only a lower limit of the DH mass, and the short-range behavior can 
be revealed. 

It is clear that both the DH structure and V q ~ q would affect the theoretical prediction of 
the DH mass. Here, we first briefly introduce the employed interaction. The Hamiltonian of 
a six-quark system in the SU(3) chiral quark model reads 



H = T + Y,^° NF + V° GE + V5 8 + V i% 

i<j 

where T denotes the kinetic energy operator of the system and V£j ONF ' , V® GE , V? s and 
represent the confinement, one-gluon exchange, pseudo-scalar chiral field induced and scalar 
chiral field induced potentials between the i-th and j-th quarks, respectively. The confinement 
potential is phenomenologically taken as 



Vg ONF = -(Xt^Ua^j + aijrfj), (2) 

which describes the long range nonperturbative QCD effect. The introduced a$ ij terms which 
take different values for different q—q pairs are called zero-point energy terms. They guarantee 
that the empirical thresholds of considered channels can more accurately be reproduced. The 
short range perturbative OGE potential is chosen to be the commonly used form [||, 10, [II, O]. 



In order to restore the important symmetry of strong interaction, the chiral symmetry, we 
introduce SU(3) chiral fields coupling to quark fields so that the medium-range NPQCD 
effects can be described || [R], |ll|]. The pseudoscalar-field-induced potentials are: 



rPS ~< ., ml 



= C(g ch , m na , A) • A(m^, A, r %3 ) (a t ■ a,) ■ {X^) f , (3) 
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and the scalar-field-induced potentials are: 



V§ = -C(g ch , m aa , A) • f 2 (m aa , A, r -) • (\?\«) f , (4) 
where the subscript / denotes that the operators in parentheses are in flavor space. The 



expressions of fi, Y, and C are shown in Ref.[ll]. To retain the important chiral symmetry 
as much as possible, we take all chiral-quark coupling constants to be the same value 

9\ _ _9_ g 2 NNn ( m i \ 2 (5) 
4vr 25 4vr \M N ) ' 1 ' 

In Eq.(3), 7r a with (a = 1, 2, 8, and 0) correspond to the pseudoscalar fields it, K, rjs and 
?7o, respectively, and rj and r( are the linear combinations of r/o and % with the mixing angle 
9. In Eq.(4), a a with (a = 1,2, ...,8, and 0) correspond to the scalar fields a', k, e and a, 
respectively. 

In solving this six-quark system problem, the first selected set of model parameters is that 



used in Ref.[ll|. This is because that with this set of parameters, almost all empirical partial 
wave phase shifts of the N — N scattering can be well re-produced, meanwhile the available 
cross sections of N — Y processes can reasonably be explained, some masses of baryon ground 
states can accurately be obtained, and some properties of baryons such as EM transition 



rates and etc. can better be understood |13|, 12]. Thus, the predicted mass of DH is based 



on a more solid ground, and the reliability of the prediction is increased. Furthermore, the 
possible range of the DH mass can be tested by shifting the values of model parameters within 
reasonable regions. 

Then, we show how to choose the model space in solving the bound state problem of a 
six-quark system. There are two types of possible configurations in studying the structure of 
DH. 

(1) Six-quark cluster configuration. 

In this configuration, the trial wave function can be expressed as the linear combination 
of differently sized basis functions: 

*(A„)y T S = Y.°i $ (AM)/ T S("i), (6) 

i 

with the basis function 

*(Am), t sfa) = <t> [ (Os) 6 , <*] x{l^ )f X c m > (7) 
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where (f> [ (Os) 6 , u>i] is the orbital wave function (oj = -^i), and xf\^ f an d X(oo) denote the 
wave functions in the flavor-spin and color spaces, respectively. This trial wave function is 
in the pure symmetry-basis-function space. In this configuration space, (Os) 6 configurations 
which describe the breath mode are considered only. 

(2) Two-cluster configuration. 

In this configuration, there exist three possible channels: AA, N3 and EE. In the frame- 
work of Resonating Group Method (RGM), the trial wave function of DH can be written 

as 

^ = a|AA) + /3|iVH)+7|EE), (8) 
with the two-cluster wave function 

| B 1 B 2 ) = A[(/)B 1 4>B 2 Xrel K-cm ]sT=0Q, (9) 

where A stands for the antisymmetrizer, 4>b 1 (2) denotes the wave function of the cluster 
B\ (2) , Xrel represents the trial wave function of the relative motion between clusters B\ and 
B2 and 1Z cm is the wave function of the total center of mass motion. This trial wave function 
is in the physics-basis-function space. The physical picture of this configuration is that in 
the compound region of two- interacting clusters (or composite particles), there might exist a 
bound state or a resonance. 

Let us define a quantity 

E H = M H - M AA , (10) 

where M# and M A a denote the mass of DH and two A's, respectively. Apparently, En < 
stands for a stable DH against weak decay. En (or Mfj) can be obtained by solving the 
Schrodinger equation in which the above mentioned potentials are employed. The results 
with different configurations are discussed in the following. 

(1) Six-quark cluster case. 

In this case, uji are taken as variational parameters. By using the variational method, one 
can minimize the Hamiltonian matrix element with respect to u>i. The resultant masses of 
DH are tabulated in Table I. 
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Table I. Eu{MeV)^ in the six-quark cluster case 



yOGE + yCONF 


311 


yOGE _j_ yCONF 

+v n + v a 


127 


yOGE + yCONF 
+ yPS + yS 


276 



{ Parameters used are those in ref. jLl[ 



It is shown that if one only employ OGE and confinement potentials, the mass of six-quark 
cluster is 311MeV heavier than those of two A's, Maa- When one additionally employs V* 
and V a , the mass of DH would decrease, but it still YllMeV heavier than Maa- However, 
when one includes all SU(3) chiral pseudoscalar and scalar fields, namely employs the other 
chiral fields in additional to the it and a fields, the corresponding mass becomes larger again, 
which is 276MeV^ heavier than Maa- This means that the couplings between a chiral field 
and quark fields cause additional attraction, which is helpful to reduce the mass of the six- 
quark system. On the contrary, the overall effect of the contributions from other SU(3) chiral 
fields provides a repulsive feature so that DH is quite hard to form. Moreover, no matter in 
which case, the mass of DH is heavier than M\\, namely DH is not a stable particle against 
to strong decay to AA and/or AfH. Therefore, a model space with (Os) 6 six-quark cluster 
structure only may not be a favored model space in studying the DH structure. 

(2) Two-cluster case. 

In RGM, to solve the bound state problem, one usually expands the unknown relative 
wave function Xrei by using locally peaked Gaussian basis functions 



Xrel 



E 



Ci Xi 



(11) 



where q 's are variational parameters | J17|| . 

Due to the Pauli principle, there exists a forbidden degree in the six-quark trial function 



m 



i 



^forbidden = A ( _L | AA) + -^= | NE) + | EE) 

^ v3 \/3 

This forbidden degree and almost forbidden degrees can be detected by examining the zero 



and almost zero eigenvalues of the normalization kernel, respectively [18|. In particular, in the 
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bound state RGM calculation, a component with the inter-cluster distance to be zero in the 
trial wave function, which is just a six-quark cluster configuration with the [6] symmetry, has 
to be included so that the behaviors of two clusters at the shorter inter-cluster distance can 
be well described and the stable and reliable solutions can be obtained. As a side-effect, the 
disturbances from the forbidden and almost forbidden degrees become serious. Sometimes, 
these disturbances would spoil the numerical calculation and produce non-physical results. 
Therefore, only after all the non-physical degrees are completely eliminated, the resultant 
energy of the bound state can be trusted. Due to the aforesaid reasons, at this moment, 
it may not be necessary to further carry out the mixing of the configurations (1) and (2). 
Moreover, in practice, eliminating the non-physical degree can be realized by performing the 
off-shell transformation. Then, carrying out the variational procedure, one can obtain the 
mass of DH or Eh- The resultant En's are tabulated in Table II. 

Table II. En(MeV) and IZ(fm) in the two-cluster case^. 





( |AA>) 


( |AA> \ 
\ + 1 NE > J 


/ |AA> \ 
+ \NE) 

V + 1 ££> / 


yOGE _|_ yCONF 


Eh 


9.41 


8.95 


6.77 


K 


1.89 


1.85 


1.68 


yOGE , yCONF 

+v n + v a 


Eh 


5.47 


-38.71 


-65.80 


K 


1.66 


0.75 


0.72 


yOGE , yCONF 

+V PS + V s 


Eh 


4.39 


4.02 


1.98 


K 


1.59 


1.57 


1.41 



f Parameters used are those in ref. jLl| and TZ denotes the root-mean-squared radius of DH. 

It is shown that in the first (V OGE + V CONF ) and third (V OGE + V CONF + V PS + V s ) cases, 
the results do not support a bound DH. In the second case (y + V + V 71 + V a ), the 
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one-channel calculation (| AA)) result (Eh = 5A7MeV) also does not support a bound state, 
but, the two- channel (| AA)+ | NE}) and three-channel (| AA)+ | NE}+ | EE)) calculations 
show a bound state with the binding energies of 38.71MeV and 65.80MeV, respectively. These 
results can roughly be explained by the interaction matrix elements. The contribution from 
OGE shows repulsive feature in the AA channel and attractive features in both NE and EE 
channels, and the net contribution from tt and a fields at the short distance demonstrates 
the more attractive character in the NE and EE channels than that in the AA channel. 
However, due to the existences of two strange quarks in the DH system, the chiral clouds 
with strangeness surrounding interacting baryons become important. Thus, in our opinion, 
all the SU(3) chiral fields should be considered. In fact, after including these fields, the 
aforesaid over-strong attraction disappears, and the results in all cases become smooth. This 
phenomenon can be understood by the matrix elements of the spin-fiavor-color operators 
of chiral fields, namely sfc coefficients, which show that the contributions from a and 7]q 
present the attractive character and those from the other mesons, i.e., n, K, r/g, a', k and e, 
are repulsive. Moreover, in the third case, although all coupled channel calculations do not 
show a bound DH, in comparison with the result in the six-quark cluster case, the inclusion 
of additional channel would reduce the mass of DH. Eventually, the resultant mass of DH is 
around the AA threshold. 

To further understand the calculated Eh, we also list the corresponding root-mean- 
squared radii (1Z) of DH in different cases in Table II. These numbers indicate that in 
the two- and three-channel calculations, if one considers OGE, n and a only, the root-mean- 
squared radii of DH are 0.75 fm and 0.72 fm for the two and three coupled channel cases, 
respectively. Therefore, DH is a bound state. In all the other cases, the resultant 1Z values 
are greater than 1.4/m, thus DH is no longer bound. 

How the values of major model parameters affect the resultant mass of DH is also studied. 
Three coupled channel results are used as samples to demonstrate these effects. When the 
width parameter b increases, the corresponding Eh value changes to a smaller value. De- 
creasing the mass of a, would make the DH mass lighter. Moreover, a smaller s quark mass 
m s corresponds to a lower DH mass, and decreasing the mixing angle 9 would just lower the 
DH mass in a very small amount. Among these parameters, the mass of a would give the 
biggest effect on Eh- 

Finally, we give the possible mass range of DH. As mentioned above, the AW scattering 
phase shifts and the NY scattering cross sections as well as the mass of DH depend on the 
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values of model parameters. Our calculation showed that except the set of model parameters 
used in Ref . [pTfl , another set of model parameters, say b = 0.53/m, m a = 600MeV, m s = 
43QMeV, 6 = -23°, A P5i a , = 987MeV, A Kj e> a , = 1381MeF, which are almost the 
values in limits, can also fit the experimental NN and NY data |19| . With this set of model 
parameters, the resultant Ejj is —6.9MeV . Further extending these parameter values to their 
physical limits with which the empirical NN and NY data cannot even been reproduced, one 
can obtain the upper and lower bounds of the DH mass. When b = 0.6/m, m a = 550MeV , 
the lower bound of Ejj is —9.1MeV. On the other direction, if one takes b = 0.48/m, m a = 
675MeV , the upper bound of DH is around 4.9MeV . 

From above calculations, one finds that in the framework of our SU(3) chiral-quark model, 
as long as one picks up a set of model parameters which satisfy the stability conditions, the 
masses of the ground states of baryons and meanwhile can be used to fit the experimental 
NN and NY scattering data, the resultant mass of DH would be rather stable and would be 
ranged in a very small region. This mass is consistent with the present experimental finding 
and reflects that the SU(3) chiral-quark model is reasonable. 

As a summery, one may have following conclusions. The six-quark system with strangeness 
being -2, J p = + (S=0, T=0) is studied in two possible model spaces. One is in a six-quark 
cluster configuration space with breath mode, and the other is in a two-cluster configuration 
space with three possible channels. It is shown that the (Os) 6 model space is not larger 
enough even the breath mode is considered. Therefore, the mass of DH in this model space 
is generally 100~300 MeV heavier than that in the two-cluster model space. The similar 
result also appears in the other six-quark system calculation, say Deltaron (d*) In the 
two-cluster configuration case, the result shows that the mass of DH is ranged from 2225MeF 
to 222>AMeV if the experimental NN and NY data should simultaneously be reproduced. It 
seems that all the SU(3) chiral fields must be considered in studying DH. In this case, the 
SU(3) chiral fields surround the baryons and make the baryons more stable and independent. 
Therefore, the interaction between two baryons becomes weaker so that it is hard to form a 
stable six-quark particle. It is also shown that the lower and upper bounds of the DH mass 
in the SU(3) chiral quark model are 2223MeF and 2237MeV, respectively. 
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